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're  to  seismic  receiver  f  >..:ict  ions  for  structure  under  PAS  (Pasadena,  CA)  are  derived  from 
az Lnut ha  1 ! \— distributed  te’.i*  seismic  I’  waves  recorded  on  Benioff  1-90  instrumentation. 

The  broad -band  *  hr ee-ci  mp.  vient  Benioff  1-90  system  is  peaked  at  one  sec  period  and 
ill. ’vs  res.  .jut  i«n  major  crustal  interfaces  from  large  Ps  conversions  seen  in  the* 
receiver  : unction  data.  !  ho  observed  bodv  wave  data  are  quite  complex  showing 
ex  eat  Iona  1 1 v  large  I's  conversions  and  scattered  waves  on  horizontal  components.  Radial 
and  tan  ;e:»t  ial  no  t ions  are  of  equal  magnitude  and  show  major  off-azimuth  converted  Ps 
waves  su.'-estin’,  large  scale  crustal  heterogeneity  beneath  the  station.'  Stochastic 
simulat  i<*;  s  of  n;  plane  lave* red  structure  show  that  geologically  unreasonable  ID  models 
arc  u  i  r."i  ■_.*  fit  the  data .  The  observed  coda  decay  yields  a  scattering  Q  estimate  of 
_■  > a  ;tC  j  S(,  ■  period  usin'.  an  energy  flux  model  for  a  propagating  plane  wave  interacting 
,  Si  at  ter  in..;  laver  over  a  homogeneous  halfspace.  Observed  and  synthetic  coda  decay 
roliows  the  theoretical  exponential  decay  predicted  by  the  model  and  is  due  entirely  to  _ 
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di; fusion  of  coda  energy  out  of  the  Layer  into  the  halfspace.  I ’AS  coda  is  compared  to 
coda  iron  deep  teleseisms  recorded  at  SCI’  (State  College,  PA)  and  it  is  seen  that  scatterin 
is  more  severe  at  PAS  as  reflected  in  higher  coda  levels  and  slower  decay  rate.  Analysis 
of  a  major  Ps  conversion  arriving  3  seconds  after  direct  P  indicates  that  a  major  crustal 
discontinuity  at  about  20  km  depth  dips  at  high  angles  to  the  north  under  the  San  Cahriel 
Mountains.  Ihis  interface  probably  represents  the  crustal  tectonic  boundary  between  tne 
ruisverso  ranges  and  the  Los  Angeles  basin. 
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?  general  objective  of  this  research  is  to  understand  the 
important  in  shaping  the  seismic  signature  of  small 
recorded  at  local  and  regional  distances.  Speciric  objec 
ire  1 )  the  identification  of  deterministic  aspects  of  the 
lid  from  small  earthquakes  and  explosions  to  allow  the 
ice  of  source  depth  and  other  source  parameters,  and  1: )  to 
:.and  the  effects  of  scattering  and  lithospheric  heteroge- 
m  the  propagation  of  high  frequency  regional  phases.  The 
id  analysis  of  deterministic"  and  "stochastic"  wave  propa- 
effect.s  is  required  to  unravel  the  complexity  of  regional 
for  the  purposes  of  event  discrimination. 


technical  Problem 


Regional  phases  from  small  events  are  affected  by  complex 
interact  ions  between  source  radiation  and  wave  propagation 
effects  due  to  structure  in  the  crust  and  upper  mantle.  Because 
observations  are  confined  to  the  high  frequency  band  {>?.  hz  )  , 
lithospheric  heterogeneity  becomes  important  in  shaping  high 
frequency  regional  phases.  Typical  wave  lengths  are  much  shorter 
than  the  total  travel  path  and  are  comparable  to  known  geologic 
~  tructure . 

An  aspect  of  the  problem  of  regional  wave  propagation  is 
examined  here  involving  the  nature  of  coda  associated  with  major 
arrivals.  Simple  wholespace  scattering  models  are  often  applied 
to  the  coda  of  regional  S  phases  to  deduce  scattering  or  anelas- 
tic  attenuation.  The  level  and  amplitude  decay  of  coda  is  a 
characteristic  of  the  data  which  seems  to  be  robust  for  particu¬ 
lar  regions  and  can  be  used  to  deduce  source  magnitude,  once  cal 
ibrated.  Factors  affecting  the  level  of  scattering  near  the 
source  and  near  the  receiver  are  also  obviously  important,  in 
yield  estimation  problems  and  waveform  modeling  for  source  par¬ 
ameters.  It  is  important,  therefore,  to  investigate  the  major 
assumptions  contained  in  these  wholespace  models  and  to  determine 
which  are  appropriate.  This  was  done  by  simplifying  the  wave 
propagation  regime  of  a  source  in  a  scattering  medium  to  that  of 
the  p>lane  wave  incident  to  heterogeneous  structure  under  a 
receiver.  Observed  coda  in  broad-band  teleseismic  receiver  func¬ 
tions  was  simulated  through  computation  of  synthetic  seismograms 
for  lb  stochastic  structures  and  through  the  construction  of 
energy  flux  models  for  wave  propagation  in  a  scattering  layer 
o'/f  r  a  homogeneous  ha  If  space.  The  results  suggest  that  the  gross 
geometry  o t  the  scatterers  is  very  important  to  the  problem  of 
coda  decay  and  coda  J  eve i  and  that  wholespace  models  signifi¬ 
cantly  misrepresent  the  coda . 


G o  ner a  L_ Methodol  ogy 

Broad  band  data  from  deep  earthquakes  recorded  at  the  sta¬ 
tion.-.  PAS  'Pasadena,  PA)  and  SCP  (State  College,  PA)  were  ana - 
lym-d  f  rom  a  point,  view  that,  the  coda  behind  direct  P  was  due  to 


c- 


wave  scattering  in  a  heterogeneous  medium  under  the  respect  i  v>; 
receivers.  The  data  at  PAS  written  by  simple,  deep 
earthquakes  show  very  complex,  high  amplitude,  and  long  duration 
coda.  An  analysis  of  Ps  conversions  in  the  first,  10  seconds  of 
the  waveforms  shows  that  structure  under  the  station  is  three 
dimensional  since  significant  particle  motions  occur  out.  of  the 
sagittal  plane  of  the  ray.  Of  interest  to  this  discussion  is 
the  nature  of  scattering  induced  by  complex  structure  and  the 
formation  of  P  coda. 

Plane  layered  stochastic  structure  models  were  constructed 
using  a  random  number  generator  with  an  exponential  spatial  cor 
relation  function  and  a  correlation  length  of  1  km  (e.  g.  ,  F'ran- 
kei  and  Clayton,  1986).  This  was  done  to  explore  the  degree  of 
ID  heterogeneity  required  to  mimic  the  observed  data.  Synthetic 
seismograms  for  an  incident  l-  plane  wave  under  the  randomly 
layered  models  were  constructed  using  the  Thompson- Haskel 1  tech¬ 
nique  . 

An  energy  flux  model  was  developed  tor  the  plane  wave 
response  of  a  scattering  layer  over  a  homogeneous  half space 
This  was  done  to  explain  aspects  of  the  ID  simulations  as  well  as 
parameterizing  the  receiver  function  data.  The  major  assumptions 
in  the  model  are  based  on  conservation  of  energy,  homogeneity  of 
the  coda  field  within  the  layer  and  leakage  of  coda  energy  into 
the  halfspace  following  a  simple  diffusion  law.  Justification 
for  these  assumptions  are  based  on  published  numerical  exper¬ 
iments  (Frankel  and  Wennerberg,  1987)  and  observations  of  coda 
behavior  in  small  regional  earthquakes  (Aki  and  Chouet,  197b; 
Dainty  and  Tokso2,  1977). 

Technical  Results 

Summarizing  the  derivation,  the  coda  amplitude  is  given  by 


=  e 

°  V  I, 


+ud/20s  (  1  -  e*“2td/0s  j  e"»t/2 


'vnere, 


Ip  =  integral  of  the  squared  velocity  of  direct  P 


td  -  h/a 


h  =  layer  thickness 
a  -  P  wave  velocity 
ex  -  circular  frequency 
0  v  scattering  0 


y  -  diffusion  constant  for  layer  boundary. 
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(1) 


Equation  (,1  i  shows  that  for  the  incident  plane  wave,  the  time 
dependence  of  the  coda  is  controlled  entirely  by  the  diffusion 
constant  of  coda  energy  diffusing  cut  of  the  layer  into  the  half- 
space .  Wholespace  energy  flux  models  have  coda  decay  controlled 
by  the  expansion  of  the  volume  behind  the  wavefront.  Single 
scattering  models  have  decay  controlled  by  scattering  0.  The 
theory  predicts  the  form  of  coda  decay  seen  in  data  and  stochas 
tic  simulations  quite  accurately.  Coda  decay  in  wholespace  mod 
els  is  only  a  resuit  of  scattering  or  anelastic  attenuation. 

Important  Find  i  rigs  .an  d  Cone  lu  s  ions 

Analysis  of  coda  in  teleseismic  receiver  functions  suggests 
that  there  are  other  mechanisms  which  control  the  formation  and 
decay  of  scattered  coda  waves  which  are  separate  from  intrinsic 
or  scattering  attenuation.  In  particular,  the  simple  geometry  of 
an  elastic  scattering  layer  over  half space  produces  coda  decay 
which  would  be  interpreted  as  an  attenuation  effect  but  is  due  to 
the  simple  redistribution  of  scattered  energy  from  the  layer  to 
the  half space.  Such  effects  can  be  studied  first  by  analyzing 
coda  from  teleseismic  events  at  a  receiver  or  array  and  then 
applying  the  parameters  of  the  layer  model  to  an  appropriate 
energy  flux  model  for  a  source  contained  within  the  layer. 

Significant  Hardware  Development 
N/A 


Special  Comments 

N/A 

Imp! i cat i ons  f o r  Further  Research 

Scattering  Q  models  are  based  on  a  number  of  assumptions  con¬ 
cerning  the  distribution  of  the  reservoir  of  energy  contained  in 
the  scattered  field.  Further  research  will  concentrate  on  com¬ 
bined  application  of  teleseismic  scattering  determin iations  and 
regional  phase  scattering.  This  combined  analysis  may  offer  con- 
staints  on  scattering  physics  not  obtainable  by  analysis  of  each 
data  set  alone.  Results  of  this  research  will  have  important 
implications  on  studies  of  regional  phase  propagation,  discrimi¬ 
nation  of  small  events,  and  yield  estimation  problems. 
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To'.eseismic  receiver  functions  for  structure  under  PAS  (P.u-aden.i.  C  \ i  ate 
e.;  ri\ed  from  azmiuthaliv  -distributed  teJeseismic  1’  waves  recorded  mi  BetiuTf  i  - 
instrumentation.  The  broad-hand  three-component  Beniofi  1-90  system  i.  peai.ed 
at  one  sec  period  and  allows  resol  u»  ion  of  major  crustal  interlaces  fmir.  large 
P>  conversions  seen  in  the  teceiver  function  data.  The  observed  bo-.h  wave  data 
are  quite  complex  showing  exceptionally  large  Ps  conversions  and  scattered  waves 
on  horizontal  components.  Radial  and  tangential  motions  are  of  equal  magnitude 
and  show  major  off-azimuth  converted  Ps  waves  suggesting  large  scale  crustal 
heterogeneity  beneath  the  station.  Stochastic  simulations  of  ID  plane  layered 
structure  show  that  geologically  unreasonable  ID  models  are  required  to  fit  the 
data.  The  observed  coda  decay  yields  a  scattering  Q  estimate  of  259  at  2  sec 
period  using  an  energy  fiux  model  for  a  propagating  plane  wave  interacting  with 
a  scattering  layer  over  a  homogeneous  halfspace.  Observed  and  synthetic  coda 
decay  follows  the  theoretical  exponential  decay  predicted  by  the  model  and  is 
due  entirely  to  diffusion  of  coda  energy  out  of  the  layer  into  the  halfspace. 

PAS  coda  is  compared  to  coda  from  deep  teleseisms  recorded  at  SCP  (State 
College,  PA)  and  it  is  seen  that  scattering  is  more  severe  at  PAS  as  reflected 
in  higher  coda  levels  and  slower  decay  rate.  Analysis  of  a  major  Ps  conversion 
arriving  3  seconds  after  direct  P  indicates  that  a  major  crustal  discontinuity  at 
about  20  km  depth  dips  at  high  angles  to  the  north  under  the  San  GabCel 
Mountains.  This  interface  probably  represents  the  crustal  tectonic  boundary 
between  the  Transverse  ranges  and  the  Los  Angeles  basin. 
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The  analysis  of  teleseismic  receiver  functions  represent'  an  inexpensive  and 


vun'enient  wax  of  imaging  major  crustal  and  upper  mantle  discontinuities  under 
i;. dated  receivers.  The  iransntisiv  ity  of  structure  under  a  three-component 
seismometer  is  inferred  from  the  timing  and  amplitude  of  IT  conversions  seen  on 
h>  rt'mtal  ground  motions  and  is  modeled  to  determine  the  location  and  velocitv 
Cs  :t  trusts  of  the  causitive  interfaces  (Burdick  and  Langston,  1977;  Langston.  1979; 
Owens.  198-1!  The  technique  has  been  particularly  useful  in  large  scale 
structure  stjdies  using  long-period  body  waves  (e.g..  Burdick  and  Langston.  1977, 
Langston  and  Isaacs,  1981;  Hebert  and  Langston.  1985)  to  determine  average 
crustal  thickness  and  is  increasingly  being  applied  to  broad-band,  high-frequency 
data  to  obtain  more  resolution  on  structure  (Owens,  1984;  Owens  et  al.  1984; 

!  9 S ”  ). 
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One  of  the  inevitable  trade-offs  in  using  higher  frequency  data  is  increased 
sensitivity  to  lateral  heterogeneity  in  crustal  structure.  In  one  sense,  this  is 
desirable  since  a  goal  of  such  studies  is  to  determine  as  much  information  about 
structure  under  the  receiver  as  possible.  However,  it  is  also  obvious  that  the 
wave  field  is  severely  spatially  aliased  through  observations  made  at  only  a 
single  surface  point.  Imaging  procedures  implicitely  rely  on  modeling  assumptions 
such  as  plane  layering  or,  at  most,  simple  curved  interfaces.  It  is  often 
observed  that  much  of  the  wavefield  is  inaccessable  to  rational  explanation 
ucing  simple  modeling  techniques  (Langston,  1979;  Owens  et  al,  1 987).  For  example, 
receiver  function  data  often  display  anomalous  wave  behavior  such  as  P-wave 
particl0  motions  which  have  significant  tangential  amplitudes. 

A  purpose  of  this  paper  is  to  examine  strategies  of  treatment  of  broad-band 
receiver  function  data  which  take  into  account  both  stochastic  and  deterministic 
scattering  effects  due  to  heterogeneous  structure.  Structure  under  the  station 
PAS  (Pasadena,  California)  will  be  the  focus  of  this  effort.  This  station  lies  in 
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'  .'.it  K  nd  u  monutn  m  <■•  n  ;i  :  >'•••:  i"d  of  linn*  n,e  «.f.>  I,  was  mitiali> 

i  ' t i \ •  t  b\  Mmr'c  ••bM,r\;i!i'-ru  t»f  complex  ;  .iriicle  m  -’ions  in  ;h<.  : -m-  r ra i :  , 

a  data  iLaucbt*  r.  1  ’  ■  lib-  data  'htc  pre\  i.  »usl>  waiuru-i  b\  l.ce«  1  ofm  i  \\  I 

attempted  to  model  the  receiwr  •  unctions  ming  ilvnamic  m\  tracing  with  tivolif 
c  >n  -.sting  !i  miogeneoii'-  laser-  separate-.*  bv  curved  30  interface".  I  t\ 
st..  ?>  w  m  part  tall  v  sv.ee-’  rsfui  in  explaining  qualitative  aspect;  the  data  tot 
eails  urtnais  tut  he  found  that  ru\  theorr  was  inadequate  to  explain  the  his. h 
amplitude  of  inferred  converted  waves  find,  the  duration  of  signal. 

In  this  paper,  the  receiver  function  data  will  be  examined  from  two  points 
\  iew.  The  first  is  from  the  standard  method  of  treating  the  data  to  infer  major 
veicctry  discontinuities  under  the  station  using  "deterministic"  aspects  of  the 
observed  data  and  simple  velocity  models.  Ps  conversions  from  telese'smic  P 
waves  are  found  to  be  unusally  large  and  are  used  to  suggest  the  existence  of 
a  large  velocity  contrast  interface  in  the  lower  crust  which  dips  to  the  north 
under  the  San  Gabriel  Mountains. 

The  other  point  of  view  is  to  ‘reat  the  data  as  resulting  from  an  unknown 
•catering  process  ami  to  attempt  to  infer  the  severity  of  wave  scattering 
under  t'-e  station.  Using  simple  measures  of  the  P  coda  amplitude  decay  along 
with  '.  U  stochastic  structure  simulations,  the  question  is  asked:  Are  the  data 
r.ms'er.t  w  ith  scattering  due  to  reasonable  plane  layered  structure?  The 
negate e  answer  for  PAS  suggests  that  such  an  analysis  can  be  generally  used 
tc  yustify  or  not  justify  a  research  effort  in  modeling  data  with  simple  plane 
1:  ‘red  structure  models.  The  severity  of  observed  scattering  under  PAS  also 
points  out  defeciencies  in  some  simple  scattering  models  and  the  need  to  develop 
appropriate  models  for  2  and  3D  stochastic  structures. 


P\S  Milieu  has  been  in  operation  since  the  mid- 1970’s  and  has  had  a  full 
mpk-mont  of  experimental  long-  and  short-period  instrumental:;,"  r'f  interest 
t  this  studs  are  the  He  muff  1-90  (Seismometer  Period  =  1  see,  Gahonometer 
Period  =  df>  see)  and,  to  a  lesser  extent,  the  Press- Ewing  30-90  systems.  The 
Her.i.d'f  1 -90  system  is  peaked  at  I  sec  period  but  records  across  a  wide 
f: e.;uer.c>  band  comparable  to  the  intermediate-period  DWWSSN  passband  (see 
Appendix )  Its  nominal  magnification  of  3000  has  allowed  routine  recording  of 
te’esci; ms  throughout  its  existence  Although  the  data  are  recorded  in  analog 
photographic  format,  the  broad  passband  of  the  instrument  potentially  allows  for 
significant  time  resolution  of  crustal  P$  conversions.  This  potential,  in 
conjunction  with  the  recognized  complexity  of  the  receiver  signal  (Langston,  1977) 
as  vs  ell  as  unusual  tectonic  problems  associated  with  the  Transverse  Ranges 
Province  of  Southern  California,  motivates  the  present  study  of  crustal 
structure  under  the  receiver. 

PAS  station  (Figure  1  )  lies  near  the  southern  boundary  of  the  Transverse 
Ranges  Province  and  the  Los  Angeles  Basin  -  Peninsular  Range  Province.  Geologic 
structure  of  the  upper  crust  is  known  to  be  quite  complex  with  major  active 
faults  separating  regions  of  diverse  rock  type.  For  example,  crystalline  rocks  of 
the  San  Gabriel  Mountains  a  few  kilometers  north  of  PAS  abut  valley  sediments 
that  attain  depths  of  up  to  10  kilometers  in  the  center  of  the  Los  Angeles  Basin 
(Yerkes  et  a!.  1965).  The  geology  of  the  Transverse  Ranges  and  the  San  Gabriel 
Mountains,  in  particular,  suggests  that  much  of  the  range  is  allochthonous  being 
thrust  over  younger  rocks  of  the  Peninsular  Range  Province. 

Hadley  and  Kanamori(  1977)  review  a  number  of  long-range  refraction  and 
travel  time  studies  for  the  area  ana  show’  that  the  Transverse  Ranges  are  a 
locus  of  both  upper  mantle  and  crustal  velocity  anomalies.  Using  the  Southern 
California  Seismic  Network,  P  delays  from  a  PKIK.P  phase  outlined  an  east-west 
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trending  A>nc  of  higher  than  average  velocitm  m  the  upper  i !v  This  /one 
"  .t'  mudied  move  e\;.:w:vei\  t  *  \  I  lumphi  ie<  et  ah  I ‘>84  t  usim  ::  :■  m:  i  rapine  im.is  i:. 

/.unique  f  •  te.o:  er  mic  i’  u  ave-  record:/  *  by  the  netw  <rk.  1  he\  deO/nrnnei 
:i,.it  the  \  n  tied  .•  -no  f  high  mantle  velocities  wu:-  contained  vv  ithir.  the 
:  r:t:v-\er.-e  ranee  province  an!  attained  depths  of  at  least  l km.  '1  !;i>  /one 

inPip:  'ted  .a-  a  region  of  man:!/  dow ni-ve  t  Dine  analogous  to  a  suKluction 
/■ne  hut  driven  In  relative  movements  of  microplates  making  up  the  San  Andrea' 
fault  system  In  this  model,  tipper  crustal  micropiate  movements  are  decoupled 
from  lower  ■c-mt  and  upper  mantle  structures  with  the  result  that  the  surface 
expression  of  the  San  Andreas  fault  is  offset  by  a  mid-crustal  horizontal  shear 
/one  from  the  fault  at  depth 

The  Transverse  Ranees  are  also  a  locus  of  change  ip.  ‘--ustnl  structure 
between  the  western  Peninsular  Ranges  -  Los  Angeles  Basin  and  the  Mojave  Block 
to  the  east  Hadley  and  Kanamori  (1977)  suggest  that  a  high-velocity  lower 
crustal  layer  comprises  about  naif  of  the  crust  in  the  Peninsular  Ranges  but 
tapers  to  oniy  a  few  kilometers  in  the  Mojave  Block.  Tomographic  study  of  Pg 
and  IT;  waves  in  Southern  California  (Hearn  and  Clayton.  19S6a:b)  support  this 
suggestion  by  showing  slower  average  velocities  in  the  Mojave  Block  relative  to 
crust  to  the  west.  Crustal  thickness  from  long-range  refraction  (Hadley  and 
Kanamcri.  19")  and  time-term  analysis  (Hearn  and  Clayton.  1986b;  suggest  that 
'C  crust  under  PAS  is  about  31  km  thick 


These  studies  demonstrate  that  structure  in  the  area  is  quite  complex.  A 
major  goal  of  the  present  work  is  to  provide  constraints  on  crustal  thickness 


under  this  important  transition  between  tectonic  provinces  as  well  as  determining 
other  structural  details.  Site-specific  information  provided  by  the  receiver 
function  technique  will  complement  these  broader  crustal  and  upper  mantle 
structure  studies. 
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Vaseform  data  from  2)  teleseism ic  earthquakes  were  obtained  from  the 
seismogram  archives  of  California  Institute  of  Technology  (Table  I)  The 
seismograms  wen.  photographed  and  enlarged  for  hand  digitization.  Waveforms 
were  digitized  at  an  irregular  sampling  interval  and  interpolated  o  an  equal 
sampling  interval  of  0.1  seconds  Processing  included  vector  rotation  of  the 
horizontal  components  into  the  theoretical  backazimuth  of  the  P-wave  arrival  to 
obtain  radial  (positive  away  from  the  source)  and  tangential  (positive  clockwise 
around  the  source)  ground  motions. 

A  source  function  equalization  procedure  was  then  performed  to  remove  the 
instrument  response  and  unknown  effective  source  function  trom  the  radial  and 
tangential  waveform  data.  In  this  procedure,  the  vertical  component  of  motion  is 
assumed  to  be  free  of  any  effect  of  near-receiver  structure  (reverberations, 
conversions)  but  contains  the  common  instrument  response  and  wave  propagation 
effects  from  the  mantle  and  near-source  region.  The  data  are  time-windowed  and 
Fourier  transformed.  The  vertical  spectrum  is  divided  into  the  horizontal 
spectra  and  then  multiplied  by  a  Gaussian  function  to  remove  high  frequency 
noise.  The  spectral  division  is  also  accompanied  by  prewhitening  the  vertical 
component  spectra  using  a  "water  level"  parameter  to  remove  spurious  spectral 
holes.  The  "water  level"  used  for  data  considered  here  was  0.1%  of  the  maximum 
of  the  vertical  component  amplitude  spectrum.  The  Gaussian  filter  used  is 
equivalent  to  a  Gaussian  pulse  in  the  time  domain  with  a  half  width  of  1  second 
(i.e..  a=1.67  in  e  c  ).  Details  and  justification  of  use  of  the  technique  can  be 
found  in  Langston(1979). 

Figure  2  shows  examples  of  waveform  data  for  events  in  the  three  major 
backazimuth  ranges.  128°.  235°,  and  315°.  Also  shown  are  PKIKP  waveforms  of  the 
1/23/76  event  which  was  also  used  in  Hadley  and  Kanamori's  (1977)  study.  This 
phase  has  unique  properties  in  a  receiver  function  study  and  will  be  used  as  a 
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er:  ul  c.  •!'.! :  amt  on  structure  under 
vN'itude  of  tangential  motions  excited  b\  the  P  wave  is  comparable  tu  radial 
i'vtioit:-.  !r,  an  ideal  rtuiiallv  stratified  earth.  I’  and  resulting  I’-.SY  conversions 
..ir  restricted  t  ■  the  sagittal  plane  containing  the  tax .  Note  also  that 

ceatia!  and  radial  motions,  or  converse!;.,  the  observed  horizontal  Component 
diffe  r  quite  substantially  in  waveform  suggesting  that  simple  instrument 
miscalibration  or  magnification  differences  cannot  give  rise  to  these  anomalous 
particle  motions.  The  Appendix  discusses  the  calibration  of  the  Beni  off  1-9  0 
system  and  shows  that  plausible  miscalibration  is  not  a  factor. 

Waveforms  for  the  PKIKP  phase  for  the  01/23/76  event  also  show  very 
anomalous  particle  motions.  This  phase  is  inmdent  below  the  crust  with  an 
incident  angle  of  about  4  degrees.  Thus,  it  is  about  as  close  to  a  vertically 
propagating  P  plane  wave  as  one  can  get  in  practice.  Horizontal  motions, 
however,  are  not  near-zero.  They  are  about  half  the  size  of  the  vertical  P 
wave  and  both  components  are  grossly  different.  Furthermore,  this  event  has  a 
backazimuth  in  which  the  observed  EW  component  is  almost  perfectly  radial,  and 
the  observed  NS  component,  tangential.  Observed  differences  cannot  be  the 
resuit  of  instrument  miscalibration  or  numerical  rotation  error.  The  data  of 
Figure  2  strongly  suggest  that  heterogeneous  three-dimensional  structure  is 
causing  large  scattering  effects  in  the  receiver  function  data. 

Because  of  the  location  of  major  teleseismic  source  zones,  the  equalized 
radial  and  tangential  component  data  were  grouped  into  three  backazimuth  groups 
and  stacked  (see  Table  1  for  groupings).  The  stacking  was  done  by  shifting  all 
traces  to  a  common  relative  time  based  on  the  P  first  arrival  and  then 
averaging  the  waveforms.  Waveforms  for  plus  and  minus  one  standard  deviation 
about  the  mean  at  each  sample  point  were  also  computed.  The  resulting  waveform 
stacks  are  shown  in  Figure  3  with  their  standard  deviations  and  a  comparison  of 
stacks  from  the  three  backazimuth  groups  are  shown  in  Figure  4.  A  display  of 


'■  within  the  \va\efunii 


thu-  t\pe  \  leiJs  information  on  the  coherency  of  '’rri-.a 
and  of  t he  level  of  processing  noise  (Owens.  1984). 

Only  the  first  15  seconds  of  the  waveforms  are  shown  in  Figure  3. 

Obviously,  the  data  of  Figure  2  show  major  arrivals  in  the  horizontal  waveforms 
for  at  least  60  seconds.  This  poses  a  dilemma  for  structure  models  that  can  be 
considered.  It  is  very  difficult  to  get  such  arrivals  from  plausible  plane 
layered  models.  This  problem  will  be  addressed  in  a  later  section.  We  will  first 
concentrate  on  major  initial  arrivals. 

Figure  3a  shows  the  waveform  stacks  for  events  from  a  backazimuth  of  128°. 
The  bounding  envelope  for  +  one  standard  deviation  is  quite  large  for  tangential 
motions  and  relatively  large  for  arrivals  after  direct  P  on  the  radial  component. 

A  phase  marked  "Ps"  on  the  radial  stack  is  observed  on  the  radial  data  from  the 
other  backazimuth  groups.  This  Ps  conversion  is  indistinct  on  the  tangential 
data  and  is  approaching  noise  level. 

The  other  two  backazimuth  groups  (Figure  3b  and  3c)  show  remarkably  large 
arrivals,  however.  The  Ps  conversion  is  roughly  half  the  size  of  direct  P  on 
the  radial  components  and  is  also  resolvably  large  on  the  tangential  components. 
Note  also  that  the  tangential  components  show  that  the  direct  P  wave  amplitude 
is  variable  within  the  noise  of  measurement.  There  is  some  hint  that  later 
arrivals  in  the  waveforms  are  coherent  but,  as  we  will  see,  are  problematical. 

Figure  4  shows  the  stacks  displayed  together  with  P  and  Ps  phases 
annotated.  The  working  hypothesis  is  that  this  secondary  phase  is  a  direct 
P-to-SV  conversion  beneath  the  station.  Its  large  size,  relative  to  direct  P,  on 
all  except  the  128°  tangential  stack  is  remarkable  and  can  be  directly  seen  in 
the  data  of  Figure  2  (Also  see  Appendix).  Figure  5  shows  particle  motion  plots 
of  the  radial  and  tangential  stacks  for  315°  and  235°  backazimuths.  Note  that 
the  direct  P  wave  conforms  to  nearly  radial  motions  as  expected  for  ideal  P 
particle  motion,  but  that  the  Ps  conversion  has  been  rotated  45°  or  more  out  of 
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the  sagittal  plane. 


The  magnitude  of  tfm  effect  i>  quite  large.  It  is  vet  \  difficult  to  produce 
.  ,ich  1’-  ar.;\,i!s  from  simple  dipping  interfaces  itial  dip  only  a  feu  degrees  or 
has e  velocitv  contrasts  typical  of  continental  crustal  models  (Langston.  !9~7; 

Lee.  19S5L  The  implication  is  that  there  is  a  major  discontinuity  under  I’AS 
which  has  high  dip  and  or  high  S  wave  velocitv  contrast. 

These  arris als  are  also  directly  evident  in  long-period  data  from  the 
station.  Figure  6  compares  observed  waveforms  for  the  1 1/29/74  event  recorded 
on  the  Benioff  1-90  and  Press- Ewing  30-90  systems  at  PAS.  The  data  have  been 
shifted  to  a  common  time  base.  The  vertical  components,  although  showing  some 
scattered  waves  in  the  coda,  are  pulse-like  and  show  one  major  P  arrival.  The  NS 
1  -90  component  is  dominated  by  the  Ps  conversion.  Direct  P  is  a  minor  initial 
arrival.  The  30-90  NS  record  shows  that  the  Ps  conversion  (arrow)  also 
dominates  the  long-period  waveform.  The  EW  components  show  the  Ps  conversion 
(arrow's)  but  direct  P,  the  first  pulse,  is  larger.  Nevertheless,  the  Ps 
conversion  is  evident  having  the  effect  of  broadening  the  initial  pulse  by  a 
factor  of  two  compared  to  the  vertical  long-period  P-wave  and  producing  a 
"shoulder"  on  this  pulse.  This  comparison  of  data  recorded  on  tw'o  different 
seismometer  systems  demonstrates  that  the  crustal  structure  responsible  for 
these  scattering  effects  is  radical  and  that  the  effect  is  not  an  artifact  of 
instrument  miscalibration. 

One  last  cor.staint  on  the  occurrence  of  this  major  Ps  conversion  can  be 
obtained  from  the  data.  Langston  (1977,  1979)  suggests  that  the  amplitude  and 
polarity  behavior  with  azimuth  of  tangential  Ps  can  constrain  the  direction  and 
magnitude  of  interface  dip.  Unfortunately,  the  phase  is  only  well  developed  in 
the  235°  and  315°  azimuth  groups  although  it  may  be  significant  that  it  is 
poorly  developed  in  the  128 0  azimuth  group.  The  PKIKP  phase  from  the  1/23/76 
event,  however,  offers  some  independent  information  in  this  regard. 


Since  tins  phase  is  neniK  vemcalh  piopagating.  an\  i’s  conversion  from  a 
dipping  interface  «  :!I  be  contained  in  the  plane  of  the  ray  and  dip  direction. 

This  is  self-evident  from  the  geometrx  (Langston.  1977 1.  Thus,  a  simple  plot  of 
particle  motion  of  the  Ps  conversion  will  lead  to  a  direct  measurement  of  the 
cop  direction.  Figure  7  displays  such  a  plot.  The  Ps  conversion  is  the  largest 
arti'al  mi  the  NS  component  and  is  easily  seen  in  the  particle  motion  plot. 

Note  that  it  is  polarized  almost  perfectly  northward.  Because  Ps  is  positive, 
relative  to  direct  radial  P  on  the  waveform  stacks  of  Figure  4.  it  is  due  to 
conversion  fron  a  higher  to  lower  velocity  as  depth  decreases.  Taking  the 
negative  polarity  of  the  PKIK.P  phase  into  account  yields  an  absolute  dip 
direction  of  northward  for  the  postulated  dipping  interface.  This  is  consistent 
with  the  tangential  Ps  polarities  displayed  by  the  235°  and  315°  stacks. 

Deterministic  Modeling  of  the  Ps  conversion 

Recent  efforts  in  modeling  receiver  function  data  have  concentrated  in 
formal  inversion  of  the  data  to  obtain  a  plane  layered  crustal  and  upper  mantle 
model  (e.g.,  Owens  et  al.  1987).  Characteristics  of  the  Pasadena  data  set  preclude 
this  approach.  Excessive  duration  and  amplitude  of  the  horizontal  component  coda 
and  large  tangential  amplitudes  all  argue  against  finding  reasonable  plane 
layered  models.  This  will  be  addressed  below.  Nevertheless,  the  distinct  Ps 
conversion  seen  in  the  data  must  have  important  implications  on  the  nature  of 
structure  under  the  station. 

Figure  8  shows  a  comparison  of  observed  and  synthetic  radial  component 
waveforms.  The  radial  stack  for  235°  azimuth  is  shown  below  a  synthetic  radial 
seismogram  computed  for  a  crustal  mode!  proposed  for  the  area  by  Hadley  and 
Kanamori(1977).  The  Thompson-Haskeli  method  (Haskell,  1962)  was  used  to  constru; 
the  synthetic.  The  crustal  model  is  shown  in  Figure  9  with  parameters  tabulated 


in  Table  2.  The  Moho  occurs  at  31  km  depth  and  produces  a  moderately  large  Ps 


lunu'isik'ii  which  arrives  4  seconds  after  direct  P  (shown  by  arrow  in  I  icure  8i. 

The  observed  Ps  conversion  is  larger  and  arrives  at  least  one  second  earlier. 

1  his  is  significantly  earlier  and  suggests  several  possibilities. 

hirst,  if  the  Ps  conversion  is  considered  to  be  from  the  Moho,  then  the 
crust  must  be  at  most  27  km  thick  if  Hadley  and  Kanamori's  velocities  are 
assumed.  Hearn  and  Clayton's  (1986a)  imaging  study  using  Pg  waves  suggested 
average  crustal  P  wave  velocities  in  the  area  of  about  6.3  km/ sec  and  a  crust al 
thickness  of  about  31  km.  It  is  possible  that  the  receiver  function  data  are 
sampling  a  local  anomaly  unresolved  by  Hearn  and  Clayton’s  data.  Alternatively,  if 
the  crust  is  31  km  thick  then  average  the  average  S  w'ave  velocity  in  the  crust 
must  be  at  least  4  km/sec.  The  average  P  w’ave  velocity  w-ould  also  be 
correspondingly  high  at  6.9  km/sec  assuming  Poisson  ratios  near  0.25,  appropriate 
for  crustal  rocks.  These  calculations  were  performed  assuming  plane  wave 
propagation  through  a  single  layer  crust  over  a  mantle  halfspace.  In  either  case, 
in  the  event  of  a  thin  crust  or  a  thicker  high-velocity  crust,  there  should  be 
an  anomalous  mass  excess  in  the  crust  and  upper  mantle  column  which  would  show 
up  in  the  gravity  field.  No  such  anomaly  is  observed  (Hadley  and  Kanamori,  1977). 

One  simple  solution  to  the  problem  is  to  accept  the  average  thickness  and 
crustal  velocities  determined  from  previous  studies  and  to  infer  a  mid-crustal 
interface  causing  the  large  Ps  conversion.  A  plane  layered  model  which  shows 
the  general  attributes  of  this  working  hypothesis  is  shown  in  Figure  9  (and  Table 
2)  and  the  corresponding  radial  synthetic  in  Figure  8.  Ps  arrivals  from  the  Moho 
are  minimized  by  making  the  structure  approximate  a  smooth  gradient  in  that 
region.  The  large  Ps  relative  amplitudes  observed  requires  a  high  S-wave 
velocity  contrast.  This,  in  turn,  implies  a  velocity  inversion  in  the  midcrust  to 
attain  the  required  large  contrast.  The  synthetic  shown  in  Figure  8  for  this 
kind  of  low  velocity  zone  (LVZ)  structure  shows  the  general  characteristics  of 
the  data  by  approximating  the  arrival  time  and  the  double-peak  character  of  the 
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The  data  may  u>  quire  that  the  1. V7.  be  '■-lie!  L\  shallow et  but  the 
gen  via!  characteristics  are  produced  b\  the  model. 

Obviously.  the  plane  levered  model  does  not  -’.-.plain  the  anomalous  particle 
motion  of  the  l’s  conversion.  A  series  of  ra\  theory  calculations  (Langston,  1977) 
..ere  performed  to  test  the  dipping  interface  model.  The  tangential  Ps  data 
require  that  the  interface  dip  northward  under  the  San  Gabriel  Mountains. 

Experience  with  such  calculations  indicates  that  there  can  be  considerable 
trade-off  between  interface  dip  and  velocity  contrast  (e.g.,  Langston,  1979). 

Waves  which  approach  a  dipping  interface  from  the  down-dip  direction  will  have  an 
effective  angle  of  incidence  which  is  larger  than  waves  approaching  a  horizontal 
interface.  This  will  produce  a  correspondiv  larger  conversion.  Of  course,  waves 
approaching  from  the  up-dip  direction  will  tend  to  have  lower  angles  of  inmdence 
with  less  conversion. 

A  number  of  dipping  interface  models  initially  based  on  the  Hadley  and 
Kanamori  crustal  model  of  Table  2  were  examined.  The  top  of  the  6.8  km/sec 
layer  was  allowed  to  dip  up  to  40°.  Two  rays  were  traced  through  the  model. 

These  were  direct  P  and  the  Ps  conversion  from  the  dipping  interface.  It  was 
quickly  seen  that,  although  it  may  be  possible  to  produce  large  Ps  conversions 
for  rays  which  approach  the  structure  from  the  down-dip  direction,  models  with 
dips  greater  than  10°  consistently  produced  low  amplitude  Ps  conversions  for 
rays  traveling  from  the  up-dip  direction.  Indeed,  for  P  velocity  contrasts  of 
6.2/ 6.8  and  5. 0/6. 8  (velocities  in  km/sec  and  assuming  a  Poisson  solid),  dips  of 
30°  resulted  in  Ps  conversions  which  had  opposite  polarities  relative  to  direct 
radial  P.  This  is  clearly  inconsistent  with  the  data  which  show  large  positive 
Ps  conversions  (Figure  4).  Thus,  interface  dip  is  required  to  be  of  the  order  of 
10°  or  less.  The  S- wave  velocity  contrast  is  also  required  to  be  greater  than  1 
km/sec.  The  calculated  Ps/P  ratio  for  up-dip  ray  incidence  is  0.19  for  the 
5. 0/6. 8  interface  at  10°  dip.  The  observed  Ps/P  ratio  for  the  315°  stack  is  0.57 


i  Figure  4'  which  represent  waxes  corning  up- dip  but  at  an  angle  from  the 
r. or ili w ar d  dip  direction.  Calculated  tangential  amplitudes  for  the  Ps  conversion 
are  comparable  to  the  radial  amplitudes  and  agree  with  the  451’  polarization 
inomalx  seen  in  the  data. 

In  summary  Ps  -•  P  arrival  times  suggest  that  a  major  discontinuity  occurs 
in  the  mid-crust  under  PAS.  Although  the  polarity  of  radial  and  tangential  Ps  is 
consistent  with  the  interface  dipping  northward  under  the  San  Gabriel  Mountains, 
dip  is  of  the  order  of  10°  or  less  and  the  S-wave  velocity  contrast  must  be 
unusually  large  (>  1  km/sec).  Qualitatively,  a  major  crustal  low  velocity  zone  can 
explain  these  observations  but  the  extremely  large  Ps/P  amplitude  ratios 
probably  imply  that  other  factors  are  affecting  the  waveform  such  as  ray 
focusing  (Lee,  19S3;  Lee  and  Langston.  19S3a;b). 


Stochastic  Structure  Modeling 

Up  to  this  point,  the  data  have  been  treated  from  a  deterministic  point  of 
view,  a  discrete  arrival  was  identified  in  the  observed  waveforms  and  modeled 
assuming  it  occurred  at  a  well  defined  interface.  The  data  show  that  the  first 
Ps  conversion  is  indeed  a  major  wave  propagation  effect  and  it  is  reasonable  to 
assume  tnat  major  early  arrivals  will  be  due  to  direct  interactions  with 
discontinuities  in  the  structure.  However,  if  the  data  of  Figure  2  are 
objectively  examined  it  becomes  clear  that  the  inferred  Ps  conversion  is  simply 
the  first  of  many  large  arrivals  in  the  horizontal  P-wave  coda.  Are  these  later 
arrivals  fundamentally  different  from  early  arrivals?  What  do  these  large 
arrivals  imply  about  the  heterogeneity  of  the  structure  under  the  station? 

The  answer  to  these  questions  are  probably  beyond  the  scope  of  this  study 
because  of  basic  limitations  in  the  data  set  and  in  knowledge  of  theoretical 
effects  of  wave  propagation  in  heterogeneous  structure.  However,  a  slightly 
different  question  can  be  asked  which  can  provide  insight  into  the  problem  of  the 
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cmia.  Can  plant-  lavered  structure  mimic  the  cuda  seen  in  the  data  and.  it' 
si-,  -a  hat  are  the  implicatum5-  ’ 

Fir  - 1  we  make  the  assumption  that  all  arriv  als  seen  after  direct  1’  in  the 
data  represent  waves  scatteied  in  structure  near  the  receiver.  This  assumption 
is  probably  poor  for  shallow  earthquake  sources  simply  because  of  known 
propagation  effects  like  near-source  surface  reflections.  Deep  events,  however, 
should  be  less  affected  by  near-source  scattering.  A  simple  geometric  argument 
can  be  made  which  requires  that  observed  scattering  occur  near  the  receiver. 

For  deep  teleseisms,  tangential  amplitudes  in  the  coda  are  seen  to  be  as 
iarge  or  larger  than  the  vertical  coda  amplitude.  Therefore,  if  these  coda  waves 
are  due  to  scattering  (P  to  P  or  S  to  P)  in  structure  near  the  source  they 
must  finally  convert  to  P  waves  to  arrive  soon  after  the  direct  P  at  the 
receiver.  If  they  are  P  waves,  they  must  have  azimuth  anomalies  greater  than 
45°  since  they  are  so  large  on  the  tangential  component.  This  contradicts  the 
near-source  scattering  hypothesis  since  the  scattered  waves  themselves  must 
come  from  teleseismic  distances  from  the  source.  The  conclusion  is  that  it  is 
possible  that  near-source  scattering  contributes  to  the  effective  teleseismic 
source  function  but  that  large  horizontal  coda  amplitudes  relative  to  horizontal 
P  must  be  due  to  scattering  near  the  receiver.  This  is  basically  the  same 
argument  used  for  the  simple  deterministic  analysis  of  receiver  functions. 

Cessaro  and  Butlerf  1 987 )  have  mads  similar  comments  in  their  analysis  of 
tangential  P  seismograms. 

Synthetic  seismograms  were  computed  for  a  series  of  receiver  models  with 
plane  layered  stochastic  structure.  The  procedure  used  by  Frankel  and  Clayton 
(1986)  was  adopted  to  generate  a  random  velocity-depth  function.  This  is 
summarized  as  follows: 

1.  A  random  series  of  normally  distributed  velocity  values,  v(z),  with  zero  mean 
and  unit  variance  was  generated  from  a  pseudo-random  number  generator.  A 
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0.25  kin  sampling,  interval  »3:  a>sapu\i  t ■  >i  a  la >ci  of  thi.i  ner 
2  The  velocity  f unction  wa>  1  "iirier  transformed  to  the  w a\ enumber  domain  u 
obtain  u'ki. 

5  An  exponential  correlation  function  Ntz)=e~‘’  •*  was  assumed  f o:  the  medium 

/\ 

xx  here  a  is  the  correlation  length.  The  wavenumber  spectrum  of  Ntz),  N(kj= 
2a  M  +  k — a “ ) .  xvas  used  to  filter  v(z). 

-b  The  filtered  velocity  spectrum  was  then  inverse  transformed  and  scaled  to  a 
wanted  velocity  variance  and  mean. 

a.  Synthetic  seismograms  were  computed  using  the  Thompson- Haskell  method 
{Haskell.  1962). 

Calculations  were  performed  with  a  Gaussian  correlation  function  as  well, 
but  the  exponential  correlation  function  proved  to  create  somewhat  more 
scattered  arrivals  since  its  spectrum  is  richer  ir.  higher  wavenumbers.  A 
correlation  length,  a,  of  one  km  was  assumed.  For  0.5  hz  waves  considered  in 
this  study,  the  corresponding  value  of  kza  is  approximately  equal  to  1.  where  k2 
is  the  vertical  component  of  the  wavenumber.  P-wave  velocity-depth  functions 
were  computed  using  the  random  number  generator  and  S  wave  velocities  derived 
try  assuming  a  Poisson  solid.  Density  was  held  constant.  Velocity  parameters 
for  the  halfspace  below  the  random  crustal  layer  were  generally  set  to  those  of 
the  lowermost  crustal  layer. 

The  interest  in  this  exercise  is  not  to  suggest  that  a  particular 
correlation  function  is  appropriate  for  earth  structure  under  PAS  but  to 
investigate  the  degree  of  scattering  required  to  mimic  the  observed  data.  The 
correlation  function  parameters  were  chosen  to  maximize  ID  scattering  effects 
within  the  chosen  frequency  band  because  the  observed  data  show  high  amplitude 
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coda.  Frankei  and  Clayton(1986)  suggest  that  a  self-similar  correlation  function 
may  be  appropriate  for  earth  structure  since  it  gives  a  frequency  independent  Q 
over  a  broad  frequency  band  and  is  consistent  with  array  data.  The  receiver 
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functu  n  data  set  I'run:  P\S  is  essenlialh  'limited  b\  the  accuiacv  of  hand 
die  dilation.  \*v  e  concentrate  on  examining  coda  formation  near  the  peak  of  the 
trument  response  u  r.icii  i«  accurately  digitized 
"I  wo  basic  \elocit>  moods  were  considered  in  the  ID  simulations.  The  first 
;  tt-isced  of  a  heterogeneous  crustal  layer  30  km  thick  with  a  mean  velocity  of 
c.s  km  >ec  and  velocity  standard  deviation  of  10‘F.  The  top  pair  of  waveforms  in 
Figure  10  ate  typical  examples  of  the  free  surface  displacements  computed  from  a 
number  of  realizations  of  the  stochastic  parameters.  Figure  9  shows  the 
corresponding  velocity -depth  functions  for  the  low  er  and  upper  pair  of 
synthetics.  The  incident  wave  time  function  assumed  was  the  time  derivative  of 
the  Gaussian  f  unction  discussed  above  in  the  source  equalization  section.  A 
quick  glance  at  there  svnthetics  and  the  data  of  Figure  2  show  that  synthetic 
coda  ie\eb  are  significantly  lower  than  those  observed  and  that  the  coda 
attenuates  quickly  with  time. 

It  might  be  expected  that  a  more  realistic  Moho  with  a  large  velocity 
contrast  would  trap  more  scattered  energy  in  the  crustal  layer.  The  middle  pair 
of  synthetics  shows  this  case  for  the  model  used  to  compute  the  upper 
syr.tr. e tics  but  assuming  a  halfspace  P  wave  velocity  of  8  km/sec.  Minor  changes 
occur  in  the  resulting  synthetics.  The  largest  change  is  to  accentuate  the  Moho 
Ps  conversion  by  about  a  factor  of  two  (arrow).  The  coda  is  largely  unaffected. 

increasing  layer  thickness  mnds  to  increase  the  duration  of  coda.  Increasing 
the  velocity  standard  deviation  to  20%  over  a  layer  60  km  thick  with  a  mean 
velocitv  of  6  km  sec  produces  synthetic  seismograms  which  start  to  mimic  the 
data  (bottom.  Figure  10;.  Large  Ps  conversions  and  reverberations  start  to 
attain  amplitudes  comparable  to  the  direct  radial  P-wave  and  coda  duration 
superficial' y  appears  to  agree  with  the  observed  data. 

Thus,  it  appea-s  that  ID  velocity  variations  in  excess  of  20%  over  a 
significant  thickness  of  the  lithosphere  are  needed  to  mimic  the  receiver  data  at 
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u  inch  xai  x  lit  vel.vitx  1  rum  a  bunt  ,v  tv  3.0  kn;  sea  li  >m  the  surface  no.,  the 
mantle  we:  a  -cale  of  about  1  '•  tv.  are  needed  ;  explain  the  data  Cleni  lx .  tins 
a  reologicalix  unreasonable. 

\n  immediate  conclusion  i<  that  other  mechanisms  of  scatter ing  are  needed 
•  e :  an  1  uhm e  body  wave  scattering  in  plane  lasers.  It  is  likely  th..:  there  is 
:  ignifican:  buds  wave  to  surface  «:ae  scattering  in  near  surface  layers  < Da i n t > 
et  a!.  I0"-!;  Aki  and  Chouet.  lc>~5;  Damtx  and  Tokscz.  !9"7).  This  is  expected 
becau-v  f  the  heterogeneous  peolopy.  I:  is  also  clear  that  2  and  3D 
structcies  are  required  in  the  area  to  produce  the  large  tangential  particle 
motions  observed  in  the  data.  Numerical  experiments  with  2D  elastic  finite 
.1  if  fere  nee  method?  show  that,  for  an  equal  amplitude  scattered  field,  velocity 
x  amnions  in  two  dimensions  are  about  half  that  needed  for  velocity  variations 
in  one  dimensionCMcLaughiin  et  ai.  1985).  Thus,  assumption  of  simple  ID  body  wave 
scattering  tor  explaining  much  of  the  receiver  function  coda  at  this  station  is 
naive. 


Energy  Flux  Models  for  Plane  W3ve  Scattering 

The  phne  layered  model  simulations  are  instructive  in  showing  deficiencies 
in  basic  modeling  assumptions  when  treating  the  receher  function  data.  However, 
would  also  be  useful  if  the  coda  data  could  be  parameterized  in  such  a  way 
:c  guide  the  analysis.  For  example,  are  there  aspects  of  the  coda  decay  which 
preclude  reasonable  plane  layered  models  in  other  situations?  Also,  what  does 
the  amplitude  of  the  coda  suggest  about  the  magnitude  of  scattering  near  a 


receiver? 

Teleseisimc  P  wave  coda  has  been  the  subject  of  a  number  of  studies  (e.g.. 

Aki,  J 9 T 3 ;  Dainty  and  Toksoz.  1977;  McLaughlin  et  al,  198.'  Levander  and  Hill.  1985; 
Frankel  and  Clayton,  1986;  Cessaro  and  Butler,  1987).  These  studies  arp  closely 
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-\  number  of  techniques  are  available  to  parameterize  the  coda  level  and 
time  dece.v  based  on  the  Born  approximation  ol  weak  single  scattering  or  of 
diffusion  of  coda  energy  for  extreme  scattering  (e.g.,  Aki  and  Chouet,  1975:  Aki, 
;wvo  bumt'.  et  al.  |9~4).  Recent  work  has  concentrated  on  simulation  studies 
usim  t  ’ r. ire  difference  acoustic  and  elastic  wave  propagation  methods(Levander 
and  Hill.  1985.  Frankel  and  Clayton,  1986)  which  implicitly  include  the  entire 
scattered  field.  Levander  and  Hill(  1985)  examined  scattering  characteristics  of  a 
rough  boundarv  between  a  surface  layer  and  underlying  halfspace  and  showed  that 
much  of  the  scattered  field  is  dominated  by  Rayleigh  wave  propagation.  Frankel 
and  Clavton  ( 1 9S6  •  and  McLaughlin  ct  a!  (1985)  examined  P-SV  propagation  in  2D 
random  media  to  examine  scattering  of  high  frequency  (f>lhz)  seismic  waves. 
Subsequently.  Frankel  and  Wennerberg(19S7)  developed  a  simple  theory  based  on 
previous  finite  difference  simulations  to  parameterize  coda  levels,  scattering 
attenuation,  and  intrinsic  attenuation  for  2  3nd  3D  scalar  wave  fields. 

The  success  of  a  number  of  receiver  function  studies  in  determining  plane 
iavered  crustal  and  upper  mantle  structure  indicates  that  the  scattered  wave 
field  mav  be  thought  of  being  composed  of  a  "coherent"  contribution  from  Ps 
conversions  and  reverberations  from  discrete  interfaces  and  a  "stochastic" 
contribution  from  smaller  scale  heterogeneities.  The  coherent  field  can  be  seen 
over  a  large  solid  angle  of  ray  paths.  The  stochastic  field  changes  quickly  with 
rev  parameter  and  rav  azimuth.  Examples  of  the  stochastic  field  are  variation 
in  tangential  P  wave  first  motions  over  the  events  of  the  235°  stack  in  Figure 
3b  as  well  as  coda  arrivals  with  long  lapse  times  from  the  first  arrival.  The 
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-uccess  of  an\  deterministic  receiver  function  stud >  depends  crmcalH  nn  the 
coherent  field  being  dominant  Howe'er,  the  incoherent  fitdd.  which  is  ii.-malis 
ignored  in  such  studies,  also  contains  statistical  information  on  the  degree  of 
betel ogeneits  in  the  structure  which  ma\  be  very  useful. 

Past  studies  have  not  entirely  addressed  energv  partitioning  and  the 
numerous  scattering  mechanisms  that  may  be  affecting  formation  of  the  full 
scattered  field  even  for  plane  wave  propagation.  This  problem  is  difficult  since 
solutions  of  the  wave  equation  for  2  and  3D  complex  structure  must  be  examined 
-a  heuristic  approach  will  be  used  here  to  develop  an  operational  theory 
appropriate  to  the  three  component  receiver  data.  The  purpose  of  doing  this  is 
to  empirically  compare  coda  levels  and  decay  between  different  receivers  for 
classification  purposes  and  to  suggest  avenues  of  research  that  will  address 
the  actual  wave  propagation  problems.  This  heuristic  approach  will  also  be  used 
to  quantify  the  differences  between  ID  structure  coda  development  and  coda 
observed  in  the  data. 

A  useful  method  of  parameterizing  the  P  coda  can  be  derived  following 
Frnnkei  and  Wennerbergf  1987).  They  examined  scalar  2D  finite  difference 
simulations  and  suggested  that  scattered  energy  behind  a  cylindrical  or  spherical 
wavefront  distributes  itself  uniformly  over  the  volume  behind  the  wavefront.  Aki 
and  Chouet  (1975)  arrived  at  the  same  conclusion  when  examining  coda  of  regional 
earthquakes.  7'his  simple  assumption  yielded  useful  formulae  for  coda  level  and 
decay  in  cases  of  strong  multiple  scattering  as  well  as  the  limiting  case  of  weak 
single  scattering.  Theory  for  an  energy  flux  model  with  plane  wave  propagation 
in  a  wholespace  will  be  developed  since  simple  wholespace  models  are  widely  used 
in  many  coda  studies.  It  is  also  useful  to  examine  a  simple  wave  propagation 
regime  to  motivate  more  realistic  scattering  models. 

Consider  first  a  source  of  plane  waves  which  radiates  two  oppositely 
propagating  plane  waves  in  a  scattering  wholespace.  The  total  instantaneous 
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energy  or  power. 


is  given  by  the  sum  of  the  direct  wave  power,  ,  and  t he 


coda  power. 


Specify  that 


(1) 


(2) 


where  t  is  time,  to  is  circular  frequency  and  Qs  is  the  quality  factor  for 
attenuation  due  to  wave  scattering.  Substitution  of  (2)  into  (1)  gives 


1  -  e 


wi/Oe 


(3) 


Coda  amplitude,  Ac,  is  related  to  the  coda  power  density,  ?£  ,  through  the 
principal  assumption  that  the  coda  energy  distributes  itself  uniformly  behind  the 
two  propagating  plane  waves.  First  we  have 

Ac =  d  ^  w 


,  d  is  a  scaling  factor.  For  P  and  S  plane  waves  d  =  i/Vp  where  p  is 
;er  .  If  cS  is  the  unit  plane  wave  area,  r  the  propagation  distance  from 
source  to  receiver  and  cr  P-wave  velocity,  then  the  coda  power  density  is 


"C  _  ~C  _  fc 

Vciurne  2  r  iS  2  a  t  <5S 


(5) 


Using  (3),  (4)  and  (5)  we  obtain 


can  be  estimated  using  the  observed  direct  wave  amplitude  and  correcting 


it  t'or  attenuation  through  the  scattering  medium.  Thus, 


-or 
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•*o 


(?) 


where  is  the  observed  direct  wave  power  and  -  r/c  .  The  factor  of  two 

comes  from  the  fact  that  two  plane  waves  are  propagating  in  the  medium  and 
contribute  to  the  scattered  field.  Plane  wave  propagation  theory  is  used  to 
obtain  the  estimate  of  direct  wave  power.  First  consider  the  integral  of  the 
square  of  the  ground  velocity.  Aft), 


A  (‘)  dt 


(8) 


Times  tj  and  t^  bound  the  direct  wave  arrival  and  are  estimated  from  the 
data.  The  direct  wave  power  is  therefore 


d  E-  =  c 


(9) 


Substitution  of  (9)  into  (7)  and  of  (7)  into  (6)  gives 


V  t 


s+w‘v/2Q, 


s  (  1  -  e 


(10) 


As  Frankel  and  Wennerberg  (1987)  show,  the  effect  of  attenuation  due  to 


scattering  determines  the  initial  level  of  the  coda  scattered  from  the  direct 


\\ but  coda  decay  with  time  E  mainly  controlled  In  the  time-dependent 


increase  in  \olunie  behind  the  wavefront  In  this  plane  wave  case,  the  1  /  V  1 
dependence  is  due  the  linearly  increasing  \olume  between  the  two  oppositely 
propagating  plane  waves. 

Relation  ( 1 )  is  useful  for  describing  the  coda  for  a  plane  wave  propagating 
in  a  thick  layer  where  coda  lapse  times  (wave  arrive  time  relative  to  the  direct 
wave)  are  less  than  the  arrival  time  of  waves  which  interact  with  the  lower 
boundary  of  the  layer  assuming  the  observation  point  is  at  the  surface.  This 
corresponds  to  lapse  times  of  less  than  5  seconds  for  receivers  on  typical 
continental  crust.  A  more  appropriate  model  for  receiver  function  coda  is 
scattering  in  a  heterogeneous  layer  overlying  a  homogeneous,  isotropic 
halfspace.  In  this  situation,  a  vertically  propagating  plane  wave  sweeps  through 
the  layer  once  on  its  wav  to  the  receiver,  reflects  from  the  free  surface  and 
sweeps  through  yet  another  time  on  its  way  back  to  the  halfspace.  Energy  is 
scattered  from  the  plane  wave  into  coda  energy. 

A  variation  of  this  problem  was  studied  by  Dainty  et  al(  1974)  and  Dainty  and 
Tok$oz(1977)  where  they  assumed  that  scattering  in  the  layer  followed  solutions 
to  the  diffusion  equation.  Also  assuming  that  all  energy  within  the  layer  was 
scattered  energy  they  obtained  the  following  analytic  solution  for  the  scattered 
energy  field  at  the  free  surface.  m(t)  (assuming  no  intrinsic  attenuation): 
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h  is  layer  thickness  and  <,  is  the  vertical  diffusivity  of  energy  through  the 
boundary  of  the  layer  into  the  halfspace.  The  coefficients.  an.  are  found  as 
solutions  of  the  following  equation: 


wluie  %  is  the  seismic  velocity  oi  the  halfspace. 

The  dominant  term  of  (11)  for  long  lapse  times  and  a  high  vertical  diffusion 
rate  can  be  shown  to  be  for  n=l.  Thus,  for  cases  where  diffusion  of  scattered 
energy  occurs  quickly,  coda  energy  decays  like 


where  a, 


This  behavior  can  be  incorporated  into  a  hybrid  model 


containing  aspects  of  plane  wave  propagation  in  a  layer  with  assumption  of 
homogeneity  of  the  scattered  field  within  the  layer. 

Consider  a  horizontal  scattering  layer  of  thickness  h  overlying  a 
homogeneous  and  isotropic  halfspace.  A  vertically  propagating  plane  wave  is 
incident  from  below,  passes  through  the  layer,  reflects  from  the  free  surface 
and  passes  back  through  the  layer  into  the  halfspace.  The  total  power  in  the 
svstem  can  be  written  as 


where  a  new  term  has  been  introduced  to  describe  the  amount  of 
instantaneous  energy  which  diffuses  out  the  bottom  of  the  layer  at  the  expense 
of  the  coda  instantaneous  energy,  .  Based  on  the  behavior  of  (11)  above,  we 
assume  that 


=  Y-r 


mmm 


which  by  simple  integration  yields 


(16) 


where  it  is  specified  that  Eo=0  at  t=0.  As  before,  the  power  in  the  direct 
wave  after  interacting  with  the  layer  is 

ED  =  ET  e' u2'd/C?s  (17) 


where  the  factor  of  two  in  the  exponent  comes  from  the  wave  passing  twice 
through  the  layer.  Obviously,  (14)  is  appropriate  for  short-duration  direct  waves 
and  times  greater  than  2tj, 

Substituting  (17)  and  (16)  into  (14)  gives 
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(18) 


Also  recognizing  the  volume  swept  by  the  plane  wave  is  now  n  6S ,  the  coda 
power  density  becomes 


(19) 


As  before,  the  total  instantaneous  energy  available  to  the  system  can  be 
estimated  from  the  observed  direct  wave  power,  En  ,  by 


(20) 


and  using  (9)  gives 
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Note  that  this  form  for  coda  amplitude  looks  superficially  the  same  as  that 
in  (10)  except  for  the  exponential  factor  of  time  in  the  numerator  of  (21). 

Indeed,  the  time  decay  of  the  coda  is  controlled  entirely  by  this  factor.  If  y  =  0 
is  assumed  so  that  no  coda  energy  can  diffuse  out  of  the  layer,  then  the  coda 
level  is  constant  for  all  time,  consistent  with  the  plane  wave  assumption  of  a 
packet  of  energy  being  homogeneously  dispersed  throughout  the  layer.  Thus,  the 
decay  of  the  coda  field  is  functionally  equivalent  to  the  leading  term  for  the 
formal  solution  (11)  and  (13)  particularly  considering  that  energy  is  proportional 
to  the  square  of  amplitude. 

Anelastic  attenuation  can  be  included  in  relations  (10)  and  (21)  as  the  factor 


e+cjtd/2QI  e-wt/20j  (22) 

where  0T  is  the  intrinsic  attenuation  of  the  medium.  The  first  exponential 
in  (22)  is  the  correction  factor  to  determine  total  energy  from  the  direct  w-ave 
and  the  second  exponential  gives  the  attenuation  of  coda  amplitude.  Note  that 
the  effect  of  coda  energy  diffusing  out  of  the  layer  given  in  (21)  is  exactly  the 
same  as  intrinsic  attenuation.  We  would  therefore  expect  that  it  would  be  very 
difficult  to  impossible  to  separate  the  two  effects  in  practice  using  the 
teleseismic  coda  data. 

An  implicit  assumption  in  developing  (10)  or  (21)  is  that  the  scattered  field 
is  of  the  same  wave  type  as  the  primary.  These  equations  are  appropriate  for, 
say,  the  scattered  pressure  field  from  an  incident  P  wave.  Even  for  simple  ID 
layered  structures,  much  of  the  scattered  field  is  composed  of  P-to-S 
conversions.  For  2D  and  3D  structures,  there  is  evidence  that  much  of  the 
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scattered  field  seen  at  the  surface  is  composed  of  low  group- velocity  surface 
waves  (Dainty  et  al  1974;  Aki  and  Chouet.  1975;  Levander  and  Hill.  1985).  Thus, 
there  is  a  procedural  problem  of  relating  observed  coda  wave  amplitude  to 
energy  since  the  wave  type  contained  in  the  coda  must  be  known  before  hand.  In 
principle,  it  is  possible  to  directly  infer  the  energy  contained  in  a  wavefield  if 
strain  observations  are  available.  However,  three-component  displacement  data 
cannot  be  used  without  assumptions  on  wavetype. 

Recognizing  these  limitations,  we  nevertheless  use  equations  (10)  and  (21)  as 
guides  to  the  analysis  of  the  thr^e-component  data.  These  equations  will  be 
useful  in  parameterizing  relative  levels  of  coda  and  coda  decay  between  isolated 
recievers  but  are  clearly  deficient  in  addressing  all  of  the  scattering 
mechanisms  which  are  probably  important  in  teleseismic  coda  development. 

Some  operational  aspects  of  examining  coda  decay  are  patterned  after 
previous  studies  (e.g.,  Richards  and  Menke,  1983;  Frankel  and  Wennerberg,  1987). 
Observed  three-component  data  for  a  single  event  are  first  narrow  band-pass 
filtered  with  a  Butterworth  recursive  filter  in  the  forward  and  backward 
directions.  The  two-pole  filter  used  had  corner  frequencies  of  0.25  and  1  hz  so 
the  following  results  are  appropriate  for  0.5  hz  waves.  Once  the  data  were 
filtered  the  intensity  of  the  direct  P  wave  at  0.5  hz  was  estimated  by  squaring 
the  signal,  choosing  tj  and  t2  (equation  8)  from  the  duration  of  large  motions  on 
the  vertical  component  and  integrating  over  this  time  interval.  The  power  of 
the  direct  wave  was  estimated  using  all  three  components  of  motion  over  the 
time  interval  inferred  from  the  vertical  component.  The  integral  of  the  squared 
velocity  used  for  equations  (10)  or  (21)  was  the  square  root  of  the  sum  of  the 
squares  of  the  integrals  from  each  component  of  ground  motion.  Each  component 
was  then  scaled  by  dividing  the  square  root  of  this  total  squared  velocity 
integral.  The  envelope  of  each  component  was  then  computed  by  forming  the 
analytic  signal  (Farnbach,  1975)  and  taking  its  modulus.  The  total  coda  time 


scrips  for  one  event  was  then  found  bs  summing  the  squares  of  the  envelopes  of 
the  three  components  at  each  time  point  and  taking  the  square  root  of  the  sum. 
Resulting  coda  envelopes  for  separate  events  were  then  averaged  to  obtain  a 
better  estimate  of  coda  level. 

Figure  1 1  shows  the  results  of  this  process  using  four  deep  earthquakes 
recorded  at  PAS  (Table  1).  Deep  events  were  chosen  to  avoid  contamination  by 
near-source  scattering  effects.  The  observed  levels  of  coda  are  very  high. 

Indeed,  an  examination  of  the  raw  data  (e.g..  Figure  2)  shows  that  much  of  the 
coda  comes  from  the  horizontal  components.  Theoretical  curves  computed  using 
equation  (10)  are  superimposed  on  the  coda  decay  curve  in  Figure  11  and  show 
that  an  apparent  scattering  Qs  of  200  to  300  is  required.  The  coda  time  decay 
appears  to  be  very  slow  and  is  roughly  consistent  with  i/Vl  found  in  this  model 
of  scrrtering. 

Figure  12  demonstrates,  however,  that  the  simple  ID  simulations  are  not 
consistent  with  coda  decay  following  equation  (10).  The  coda  curve  for  the  "20%" 
model  was  constructed  by  stacking  10  vertical  and  horizontal  component 
realizations  (20  time  series  in  all)  of  models  which  had  a  velocity  standard 
deviation  of  20%  and  a  layer  thickness  of  60  km.  The  *'10%"  curve  was  obtained 
by  stacking  9  vertical  and  horizontal  component  realizations  (18  time  series)  for 
models  which  had  a  velocity  standard  deviation  of  10%  and  a  layer  thickness  of 
30  km.  The  observed  coda  decay  is  linear  on  the  logarithm  plot  and  falls  off 
much  faster  than  implied  by  (10).  The  linear  fall-off  is  consistent  with  the 
scattering  layer-over-halfspace  model  where  coda  energy  diffuses  out  of  the 
layer  into  the  halfspace  governed  by  equation  (15).  The  ID  simulations  included  no 
effect  of  anelastic  attenuation. 

Figure  12  also  shows  least-squares  linear  fits  to  the  coda  decay  to  obtain 
Qs  and  Y  in  equation  (21).  The  slope  of  the  log-coda  curve  yields  Y  and  the 
zero  time  intercept  can  be  used  to  directly  solve  for  Qs.  The  standard 


de%  u.tion  of  the  least-squares  fit  was  also  used  to  estimate  allowable  Qs 
\.,n.;ti'.'n  by.  adding  and  subtracting  the  standard  deviation  from  the  zero 
intercept  time  to  find  a  lower  and  upper  bound  of  Qs.  respectively.  These 
values  are  displayed  Table  3.  Coda  from  the  simulations  show  that  the  diffusing 
layer  model  correctly  predicts  the  form  of  coda  decay  although  the  model  does 
n't  'OMnallv  treat  the  scattering  mechanism  of  P  to  S  conversions  within  the 
1  ,ver  The  decay  rate  is  very  sensitive  to  the  velocity  standard  deviation  but 
O.  estimates  are  surprisingly  the  same,  within  the  error  of  determination. 

It  is  interesting  to  compare  resuhs  for  PAS  with  those  from  another 
station  to  set  an  appreciation  for  the  level  of  scattering  implied  by  the  data. 

Three  deep  events  recorded  on  the  broad-band  DW\\SSN  system  at  SCP  (State 
College.  PA)  were  analyzed  in  the  same  way.  Event  parameters  can  be  found  in 
Table  4  and  the  data  are  displayed  in  Figure  13.  The  Benioff  1-90  and 
intermediate-period  DWWSSN  systems  are  sufficiently  similar  for  the  purposes  of 
this  comparison,  particularly  since  the  same  band-pass  filter  was  used  on  the 
data. 

Figure  14  compares  the  coda  decay  curves  for  PAS  and  SCP.  Structure  under 
SCP  is  seen  to  be  simplier  than  that  at  PAS  (Langston  and  Isaacs,  1981;  Ammon, 
19SS  personal  communication)  and  gives  rise  to  low:er  amplitude  Ps  conversions  as 
well  as  coda.  Coda  decay  for  SCP  is  twice  as  fast  than  that  observed  for  PAS 
(Table  3).  Qs  is  found  to  be  lower  for  PAS  with  use  of  equation  (21)  giving  a 
value  of  239  compared  to  582  for  SCP. 

Discussion 

The  scatterine  laver-over-halfspace  model  reproduces  the  principal  behavior 
of  the  ID  structure  simulations  (Figure  12)  and  is  consistent  with  coda  decay  in 
the  PAS  and  SCP  data.  The  simple  assumptions  of  homogeneity  of  the  coda  field 
and  diffusion  of  energy  into  the  halfspace  seem  to  describe  the  basic  mechanisms 
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of  .uJ.i  formation  and  is  consistent  with  previous  observations  of  the  behavior 
of  data  and  2D  simulation  studies 

\ki  and  Chouetl  1 975 !  estimated  the  diffttsivity  of  the  lithosphere  in  Japan 
and  California  using  a  diffusion  model  of  coda  formation  applied  tu  local 
earthquake  data.  They  found  high  diffusion  rates  having  the  effect  of 
homogenizing  the  coda  field  behind  the  wavefront.  Frankel  and  Wennerbergl  19S7 ) 
took  these  ideas  further  by  examining  the  coda  field  in  finite  difference 
simulations  and  constructing  a  simple  energy  flux  theory  to  explain  the  formation 
coda.  Although  the  assumptions  of  homogeneous  coda  and  diffusive  energy  flow 
across  the  layer  boundary  are  reasonable,  the  actual  mechanisms  of  coda 
formation  are  not  directly  addressed  in  an  equation  like  (21)  which  leads  to  the 
problem  of  estimating  coda  energy  from  an  unknown  wavefield. 


Much  of  the  coda  in  the  ID  simulations  is  a  product  of  P  to  S  conversions 
and  reverberations.  The  energy  scattered  into  S  waves  is  obviously  a  function 
of  ray  parameter.  As  the  direct  wave  incidence  angle  increases,  more  P  to  S 


conversions  will  occur.  This  can  be  verified  directly  by  calculation  but  can  be 
easily  seen  in  the  behaviour  of  the  conversion  coefficient  at  a  boundary.  Thus, 
it  can  be  expected  that  coda  fall-off  and  levels  will  change  for  waves  of 
different  incidence  angle  if  ID  structure  is  appropriate.  P  to  S  scattering  in 
2D  structures  is  more  complex  (Frankel  and  Clayton,  1986;  McLaughlin  et  al  1985) 
but  appears  to  become  less  sensitive  to  incidence  angle.  P  and  S  to  Rayleigh 
scattering  is  probably  a  major  component  of  the  coda  field  at  relatively  low- 
frequencies  (dnz)  (Aki  and  Chouet,  1975).  These  scattering  mechanisms  may  control 
the  coda  formation  in  the  data  presented  here.  In  terms  of  the  application  of 
equation  (21)  the  problem  amounts  to  estimating  d,  the  scaling  factor  relating 
energy  density  to  wave  amplitude  in  equation  (4).  Even  for  ID  structure  d 
depends  on  incidence  angle  and  includes  the  free  surface  receiver  functions  (e.g., 


Helmberger,  1968).  It  is  of  some  interest  to  examine  the  energy  partitioning  in 
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the  ;uii.t  of  the  simulations  and.  making  some  simple  assumptions,  the  partitioning 
seen  in  the  data. 

As  an  approximation,  consider  the  coda  power  being  composed  of  S-wave,  El¬ 
and  P-wave,  E..-  ,  powers 


-  r  x  c 


Also  define  the  energy  partitioning  coefficient  by 


For  plane  wave  propagation 


=  c  a  j 


where  I  and  I  are  the  estimated  integrals  of  squared  velocity  for 

s  p 

S-wave  and  P-wave  motions,  respectively.  The  S-wave  velocity  is  given  by  /5  .  For 
ID  structure  models  and  for  incident  P  waves  of  small  incidence  angle,  S-waves 
occur  primarily  on  the  horizontal  component  and  P  on  the  vertical.  The 
respective  wave  integrals  can  therefore  be  directly  estimated  using  equation  (8) 
by  performing  the  integration  over  the  filtered  and  squared  waveform  from  the 
end  of  the  direct  wave  arrival  to  some  reference  time  in  the  coda.  This  was 
done  for  the  waveforms  obtained  from  the  10%  and  20%  ID  simulation  models.  Both 
mode’s  give  similar  results  where,  for  the  ray  parameter  considered  (0.06  sec/km) 
in  constructing  the  synthetics,  roughly  70%  of  the  scattered  energy  occurs  as 
S-wave  energy  (  A  value  of  0.7  +  0.5  was  obtained  for  both  simulations  using  the 


individual  waveforms  of  each  model  realization.)  Changing  the  ray  parameter  to 
0.04  sec/km,  appropriate  for  source  distances  of  85°  reduces  the  S-wave  energy 
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to  50',  and  less  tor  the  synthetics.  The  free  suiface  eft  eel  is  assumed  tu  be 

the  same  for  both  P  and  S  for  these  loss  angles  of  incidence. 

Small  differences  in  correcting  the  data  for  wave  propagation  effects  are  of 

little  consequence  to  this  discussion  since  wave  t\ pes  in  the  observed  coda 

data  are  largely  unknown.  We  treat  the  observed  data  in  the  san  ?  way  where 

the  S-  and  P-wave  integrals  are  defined  as: 

i  I  z 

:  =  v s 

S  -  I — v  , 


where  the  subscripts  Z,  N.  and  E  denote  the  component  of  ground  motion. 
Assuming  only  S  and  P-wave  partitioning  in  the  PAS  data  yields  a  partitioning 
coefficient  of  1.7  +  0.4,  a  factor  of  2  to  3  greater  than  expected  compared  to 
the  ID  simulations.  This  result  is  consistent  with  the  coda  being  comprised  of 
low  group-velocity  surface  waves  scattered  from  incident  P  and  S  waves. 
Instantaneous  energy  will  be  proportional  to  the  group  velocity  so  that  assuming 
a  higher  S  velocity  in  (26)  will  cause  the  energy  to  be  overestimated.  These 
observations  are  consistent  with  observations  of  the  coda  at  arrays  (Aki,  1973; 

Aki  and  Chouet,  1975)  and  from  theoretical  wave  propagation  calculations  (e.g., 
Levander  and  Hill,  1985).  Powell  and  Meltaer  (1984)  also  found  direct  evidence  for 
a  high  level  of  scattering  under  southern  California  in  their  study  of  coherency 
across  the  SCARLET  array. 

It  thus  appears  that  coda  level  and  coda  decay  at  PAS  is  inconsistent  with 
plausible  ID  earth  models.  The  observed  data  show  slow  coda  decay  implying  a 
relatively  long  dwell  time  of  coda  energy  in  the  crust  as  well  as  high 


amplitudes.  The  high  amplitudes  are  consistent  with  the  coda  field  being 
primarily  composed  of  scatrered  surface  waves.  Even  the  coda  decay  seen  at 
SCP  implies  unreasonable  ID  structure  since  the  data  imply  virtually  the  same 
attenuation  and  decay  as  the  20%  velocity  model  (Table  3). 
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1  hfsc  aspects  of  the  receiver  function  data  can  be  routinelv  quantified  in 
■  the:  data  set-  to  its  tivau  an  interpretation  of  structure  under  a  receiver 


I:  the  tecei'e:  t'uncti  n  c>via  data  show  tendencies  that  are  inconsistent  with 
mm;  10  models  then  inversion  of  selected  phases  at  long  lapse  times  (>10  sec) 
an  the  first  arrival  or  inversion  of  the  entire  waveform  becomes  suspect.  This 
>c'i:h  :  not  'urprising  when  one  simply  looks  at  the  anomalies  contained  within 
the  data  but  the  theoretical  treatment  presented  here  can  help  quantify  both 
the  itrew  characteristics  of  the  observed  data  and  the  justification  of  a 
particula:  modeling  sirateev 

The  large  size  of  the  initial  inferred  Ps  conversion  with  anomalous  particle 
motion  implies  the  existence  of  an  interface  which  dips  to  the  north  under  the 
Sar.  Gabriel  Mountains  with  a  gentle  dip  of  10°  or  less.  The  large  relative 
amplitude,  seen  also  on  long-period  data,  implies  an  S-wave  velocity  contrast 
greater  than  1  km/sec.  It  is  conceivable  that  focussing  due  to  a  curved 
interface  accentuates  the  Ps  conversion  relative  to  direct  P.  However,  it  is 
difficult  to  produce  such  effects  unless  the  interface  radius  of  curvature  is 
comparable  to  the  depth  of  the  interface  (Lee,  1983)  so  that  a  caustic  is  formed 
at  the  surface. 

The  northward  dip  of  the  interface  suggests  that  it  is  a  major  structure 
associated  with  the  southward  overthrusting  of  the  San  Gabriel  Mountains.  If 
tne  amplitude  of  the  conversion  is  due  to  large  velocity  contrasts,  then  a  low 


v  elect: »  torn.?  L  required  at  mid  to  deep  crustal  levels.  It  is  interesting  to 
note  that  this  low  velocity  zone  occurs  just  under  the  seismogenic  zone  of  the 
region  and  may  be  the  seismic  signature  of  the  decoupling  zone  of  upper  crustal 
and  lower  crustal-upper  mantle  microplates  (Humphries,  1984). 


The  receiver  function  data  set  for  PAS  suggests  the  scattering  is  occuring 
in  a  highly  heterogeneous  crust.  Broad-band  Benioff  1-90  data  from  teleseisnis 
show  anomalous  tangential  particle  motions  and  a  high  amplitude  coda  which 
decays  slowly.  Initial  portions  of  the  radial  and  tangential  receiver  functions 
show  a  coherent  inferred  Ps  conversion  which  displays  a  polarization  anomaly  of 
45°  for  most  data.  Using  the  amplitude,  polarity,  and  timing  of  this  phase  seen 
in  stacks  of  the  data  and  from  a  direct  observation  in  an  incident  PK.1K.P  phase, 
a  high  S- wave  velocity  contrast  (>1  km/sec)  interface  is  inferred  at  approximately 
20  km  depth.  The  interface  dips  less  that  10°  to  the  north  and  appears  to  be  a 
major  structure  associated  with  southward  overthrusting  of  the  San  Gabriel 
Mountains. 

Observed  coda  level  and  decay  was  examined  using  two  methods.  One  was 
direct  simulation  of  ID  stochastic  structures.  Plane  wave  synthetic  seismograms 
were  computed  for  random  plane  layered  models  with  an  exponential  correlation 
function  and  with  10%  and  20%  standard  deviations  in  velocity.  The  PAS  data 
showed  larger  scattering  effects  than  the  simulations  indicating  that  geologically 
unreasonable  ID  models  are  required  to  explain  the  coda  data.  The  ID  models  also 
are  obviously  deficient  in  explaining  the  degree  of  off-azimuth  scattering  seen 
in  the  data. 

The  other  method  consisted  of  examining  coda  behavior  using  an  energy  flux 
model  developed  for  a  scalar  plane  wave  incident  on  a  scattering  layer  over  a 
homogeneous  halfspace.  A  scattering  layer  model  was  considered  since  it  is  likely 
that  major  velocity  perturbations  are  largely  confined  to  the  crust.  Two 
fundamental  assumptions  were  made  to  develop  the  model  and  were  based  on 
previous  empirical  observations  of  the  behavior  of  earthquake  coda  and  numerical 
experiments.  It  was  assumed  that  the  coda  field  distributes  itself  homoge¬ 
neously  within  the  layer  and  that  coda  energy  diffuses  across  the  layer- 


halfspace  boundary  Coda  decay  is  seen  to  be  controlled  entirely  by  the 
diffusion  constant  of  energy  flow  across  the  layer  boundary.  Synthetic 
seismograms  from  the  ID  simulations  show  that  the  simple  energy  flux  model 
explains  the  form  of  coda  decay.  One  implication  of  this  model  is  that  the 
diffusion  effect  is  indistinguishable  from  anelastic  attenuation.  Thus,  it  is 
likelv  that  teleseismic  coda  data  cannot  be  used  to  estimate  local  anelastic 
attenuation. 

PAS  and  SCP  data  from  selected  deep  earthquakes  were  analyzed  using  this 
model  and  it  was  found  that  PAS  had  a  lower  scattering  Qs  (-239)  compared  to 
SCP  (-582)  and  that  the  coda  decay  for  SCP  was  twice  as  fast  as  that  for  PAS. 
The  absolute  values  of  scattering  Qs  obtained  with  the  model  are  subject  to 
assumptions  on  the  types  of  waves  contained  with  the  wavefield  and  probably 
represent  lower  bounds  to  the  actual  Q  values.  The  comparison  between  the  two 
stations  shows  that  scattering  is  lower  in  a  tectonically  quiescent  area  with 
less  variable  geology  as  expected. 

An  analysis  of  energy  in  the  horizontal  and  vertical  components  of  the  ID 
synthetics  and  the  PAS  data  suggests  that  much  of  the  energy  contained  in  the 
observed  coda  is  from  scattered  surface  waves. 

Finally,  the  observational  techniques  proposed  here  can  be  used  to  justify 
or  not  justify  the  use  of  ID  plane  layered  inversion  models  or  the  analysis  of 
isolated  phases  late  in  the  coda  of  receiver  function  data.  The  PAS  data  are 
clearly  inconsistent  with  geologically  plausible  ID  models  as  are  the  SCP  data. 
Further  research  is  needed  to  quantify  the  nature  of  2  and  3D  scattering  in 
causing  teleseismic  coda  and  in  quantifying  the  application  of  simple  energy  flux 
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Nominal  instrument  constants  for  the  Benioff  1-90  system  are  pendulum 
period  of  1  sec.  gahenumeter  period  of  90  sec.  damping  constants  of  1  for  both 
the  seismometer  and  gal  mieter.  a  coupling  constant  of  0.05,  and  magnification 
of  5000  (H.  Kanamori.  personal  communication  1987).  A  calibration  of  the  system 
was  started  in  1962  but  was  never  totally  completed.  Calibration  of  the  vertical 
component  showed  a  peak  magnification  of  2700,  10%  under  nominal  specifications. 
Experience  with  the  system  suggests  that  instrument  constants  are  good  to 
about  50%.  Because  calibration  of  the  instruments  can  affect  the  results  of 
rotation  of  the  data  and  the  source  equalization,  it  is  of  some  interest  to 
examine  the  results  of  errors  in  the  instrument  constants.  It  will  be  assumed 
that  the  receiver  response  is  ideal  consisting  of  motions  confined  to  the 
sagittal  plane  containing  the  ray.  Vector  rotation  of  the  horizontal  displacement 
components  to  obtain  radial  and  tangential  ground  motions  follows 


(9)  -  u-d)  sin[9) 


-ClW 


where  the  subscripts  R  and  T  denote  radial  and  tangential  motions. 


respectively,  and  the  subscripts  N  and  E  denote  north-south  and  east-west 
motions,  respectively.  5  is  the  backazimuth  angle  to  the  source  from  the 
receiver.  If  the  respective  instrument  responses  are  given  by  in(t)  and  ie(t), 
and  R » t )  is  a  common  radial  response  for  plane  layered  structure,  then 

9A  (u  =  Rdi  *  {  ijvt  cos2e  -  id)  s:n2s }  ) 
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responses  then  the  radial  motions  will  be  little  affected  since  the  net  response 
will  the  average  of  the  two.  If  magnifaction  is  known  to  within  some  factor  6 
such  that 


(tj  =  (  1  +  e  )  i(l) 


:e(U  =  (  i  +<  )  i(t)  I 

where  i(t)  is  the  nominal  response  then,  assuming  a  backazimuth  of  45°, 

1 1  =  P  I*'.  *  ;  or  /  i  j.  ,  \  '\ 

«'i'l  » \<-l  (,  i  t  i 


(A4) 


uTa;  =  Rd)  *  iw  { i  e 


(A5) 


Thus,  the  tangential  component,  in  the  worst  case,  will  approximately  be  30% 
of  the  radial  component  if  magnification  is  only  known  to  30%  and  will  look 
identical  to  the  radial  component  in  waveshape.  The  data  for  PAS  (Figure  2)  snow 
extreme  differences  in  waveshape  between  the  horizontal  components  which  cannot 
be  due  to  magnification  errors. 

The  equations  for  electromagnetic  seismographs  (Hagiwara,  195S)  were  used  to 
estimate  the  difference  in  instrument  responses  if  30%  variations  in  the 
instrument  constants  are  assumed.  Figure  A1  displays  the  results  for  amplitude 
spectra.  30%  changes  in  the  galvonometer  period  and  damping,  and  the  seismometer 
period  and  damping  were  assumed  relative  to  the  nominal  response.  Theoretical 
responses  were  calculated  in  the  frequency  domain  and  inverse  Fourier 
transformed  to  obtain  impulse  responses.  Using  equation  (A3)  as  a  guide,  the 
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perturbed  responses  were  then  subtracted  front  the  nominal  response  and  then 
Fourier  transformed  to  obtain  the  amplitude  spectrum.  Thus,  the  four  curses 
below  the  nominal  1-90  response  seen  in  Figure  A 1  are  the  amplitude  spectra  of 
the  differenced  impulse  responses.  They  can  be  considered  numerical  derivatives 
of  the  instrument  response  if  divided  by  0.3.  According  to  equation  (A2).  these 
curves  would  also  be  the  spectra  of  the  tangential  impulse  response. 

A  change  in  the  galvanometer  period  or  damping  results  in  a  response  2  to  3 
orders  of  magnitude  lower  than  the  nominal  response  in  the  band  centered  about 
1  hz  (Figure  Al).  Thus,  it  is  not  likely  that  errors  in  these  parameters  will  be 
of  any  consequence  in  the  data.  Changes  of  30%  in  seismometer  period  produces 
a  tangential  impulse  response  about  a  factor  of  4  lower  than  the  nominal 
response  and  looks  nearly  identical  to  the  nominal  response.  A  change  in 
seismometer  period  appears  as  a  change  in  magnification.  The  tangential 
waveform  would  differ  by  only  a  constant  compared  to  the  radial  waveform.  A 
change  in  seismometer  damping,  however,  has  the  greatest  change  in  the  shape  of 
the  spectrum.  A  90°  phase  shift  is  evident  at  1  hz  which,  for  band  limited  data, 
would  make  the  tangential  motion  appear  Hilbert  transformed  compared  to  radial 
motion.  Fortunately,  this  response  is  1  to  2  orders  of  magnitude  lower  than  the 
nominal  response.  In  summary,  plausible  changes  in  the  instrument  constants  for 
the  1-90  system  cannot  explain  the  anomalous  particle  motions  seen  in  the  data. 

It  cannot  be  discounted  that  magnification  errors  may  occur  in  the  data.  An 
empirical  test  w-as  made  by  comparing  the  ratio  of  north-south,  east-west,  and 
vertical  amplitudes  of  the  first  P  pulse  observed  in  the  1-90  data  with  that 
seen  in  the  30-90  data.  P-wave  data  for  the  12/28/73,  03/23/74,  and  11/29/74 
events  were  used.  Considering  that  the  pass-bands  of  the  two  instruments  are 
different,  amplitude  ratios  of  the  different  ground  motion  components  between 
instruments  were  within  20%  of  each  other. 

Finally,  the  data  can  be  used  in  a  the  test  proposed  by  Langston  (1979)  to 


dem  nsiiate  that  major  ot'l -azimuth  arris  uis  occur  on  both  hwii/.untai  component' 
u-i  events  with  different  hackazimuths.  Figure  A2  shows  the  lira  10  seconds  o! 
die  I’  waveforms  for  the  02  01  73  and  12  28  73  events  (also  see  figure  2). 

!\.!a:  ities  base  been  adjusted  to  make  the  waveform  comparison  clearer  and  the 
I’-’.v.T-es  have  been  aligned  in  time.  The  arrow  in  the  middle  pair  ol  plots  shows 
the  location  of  the  large  Ps  conversion  studied  in  the  main  body  of  the  paper. 

1:  clearly  occurs  on  the  north-south  component  for  the  02/0 1  '73  event  and  on 
the  east-west  component  of  the  12,  28 '73  event.  Likewise,  the  corresponding 
east- west  and  north-south  components  (right  pair  of  waveforms)  show  similar 
waveforms  between  events  without  the  major  arrival.  This  comparison  shows  that 
both  horizontal  instruments  behaved  in  a  similar  fashion  for  the  same  wave 
propagation  effect. 
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Table  2  -  Plane  Layered  Crustal  Modeb 
Laser  V  n  (km  sec)  Vs  (km  sec)  Density  (env'cc)  Thickness  (km 
Hadlev  and  Kanamori  Model 
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4.27 

2.95 

8 

7.6 

4.39 

3.0 

9 

8.0 

4.62 

3.2 

10 

8.3 

4.8 

3.35 

Table  3  -  O.  and  1  Determinations  for  Coda  Stock? 
Least-Squares  Fit 

C .via  Stack _ f n te rceot  Slope  Stand. Dev.  Q,. _ >  High  Q  Low  O 


Simulation 

-0.070 

-0.050 

0.079 

544 

0.229 

784 

380 

Simulation 

-0.985 

-0.009S 

0.057 

584 

0.045 

761 

450 

•\S  L>ata 

-0C89 

-0.0047 

0.075 

239 

0.022 

337 

169 

CP  Data 

-0.984 

-0.0092 

0.105 

582 

0.043 

944 

359 

'Scattering  Q  determinations  where  found  using  equation  (21)  for  a  30  km  thick 
scattering  layer  with  an  average  P  velocity  of  6  km/sec.) 
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Table  4  -  Event  Parameter?  for  State  College  Data 


Origin  Lat.  Lon.  Mb  Depth  Dist.  BAZ 


12/21/83  12:14:18 
08/08/85  16:35:58 
08/12/85  04:36:43 


63.2W 
113.4E 
!  1 7.2E 


602 
5.7  596 
5.6  599 


70.0  166 

144.0  341 

143.8  334 


AM 


i  mure  Captions 

Figure  1:  Sketch  map  of  southern  California  showing  the  location  of  PAS  station 


and  nta.'or  faults  of  the  area. 

Figure  2:  Selected  data  from  deep  earthquakes  (see  Table  1)  recorded  on  Denioff 
1-90  instrumentation  at  PAS.  2.  N'S,  and  EW  waveforms  denote  observed 
vertical,  north-south,  and  east-west  components.  R  and  T  waveforms  are  the 
result  of  rotating  the  observed  horizontal  components  into  the  theoretical 
backazimuth  of  the  ray.  The  distance  and  backazimuth  angle  are  given,  in 
order,  in  the  parantheses  to  the  right  :  each  event's  date.  Note  the  time 
scale  difference  for  event  01/23/76. 

Figure  3:  Stacks  of  source  equalized  radial  (top)  and  tangential  (bottom)  data  for 
the  three  azimuthal  groups  considered  in  the  text.  The  average  and  +  one 
standard  deviation  waveforms  are  shown  in  each  panel.  P  and  Ps  arrivals 
are  annotated. 

(A)  128  Degree  group 

(B)  235  Degree  group 

(C)  315  Degree  group 

Figure  4:  Comparison  of  the  stacked  radial  and  tangential  equalized  waveforms 
for  the  three  backazimuth  groups.  Note  the  azimuthal  dependence  of  the  Ps 
conversion  on  the  radial  components.  The  waveforms  have  been  shifted  in 
baseline  for  viewing  purposes. 

Figure  5:  Particle  motion  plots  for  the  315  and  235  degree  equalized  waveforms. 

The  waveform  data  are  displayed  above  the  radial  (R  or  RADI)  -  Tangential  (T 
or  TANG  I  particle  motions.  Arrows  are  shown  on  the  particle  motion  plots 
every  C.5  sec.  Waveform  data  included  within  the  brackets  are  plotted 
below.  Maximum  amplitude  for  each  waveform  is  shown  to  right  of  the 
waveform.  Note  that  the  P  waves  in  both  cases  are  polarized  in  the 
expected  ray  direction  but  the  inferred  Ps  conversions  have  a  polarization 
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annnv.ilv  ot  4a  decrees. 

Figure  6:  Comparison  of  data  from  the  11/29/74  event  recorded  on  1-9Q  (bottom', 
and  30-90  instrumentation  at  PAS.  The  arrows  show  the  location  of  the  re¬ 
conversion  discussed  in  the  text  for  both  data  sets.  Note  particularly  the 
extreme  amplitude  of  the  Ps  conversion  relative  to  the  first  P  arrival  on 
both  IPS  and  LPN  components. 

Figure  7;  Particle  motion  plot  of  the  horizontal  data  for  the  01/23/76  event. 

Data  included  within  the  bracket  are  plotted  below  with  arrows  occurring 
every  0.5  seconds.  The  Ps  conversion  is  polarized  almost  perfectly 
northward  and  indicates  the  direction  of  dip  of  the  causitive  interface. 

Figure  8:  Comparison  of  synthetic  radial  waveforms  for  the  two  models  of  Taole 
2  with  the  235  degree  radial  waveform  stack.  The  arrow  for  the  H-K 
(Hadley-Kanamori  model)  waveform  shows  the  location  of  the  Ps  conversion 
from  the  Moho.  The  arrow  for  the  LVZ  (Low  Velocity  Zone  model)  waveform 
show's  the  location  of  the  Ps  conversion  produced  at  the  base  of  the 
crustal  low  velocity  zone.  The  baselines  of  the  synthetic  waveforms  have 
been  shifted  for  viewing  purposes. 

Figure  9:  Velocity  depth  functions  for  the  Hadley-Kanamori  model  and  the  LVZ 
model  used  in  constructing  the  synthetic  seismograms  of  Figure  8.  Also 
shown  are  the  stochastic  ID  models  used  in  the  calculation  of  the  synthetics 
displayed  in  Figure  10.  Note  that  velocity  and  depth  scales  are  different 
among  the  plots. 

Figure  10:  Typical  plane  wave  synthetic  seismograms  computed  for  three 
realizations  of  ID  stochastic  velocity  structure.  Vertical  and  radial 
displacement  components  are  shown.  The  top  traces  are  for  a  model 
containing  a  layer  30  km  thick  with  an  average  velocity  of  5.5  km/sec  and  a 
standard  deviation  10%  of  the  average.  The  center  pair  of  synthetics  are 
for  the  same  crustal  model  as  the  top  pair  but  with  the  addition  of  a  high 


\  e i o c i t \  (8.0  km  sec)  halfspace.  The  arrow  points  to  the  Molio  Ps  conv  ersion 
on  the  radial  component.  The  lower  pair  of  synthetics  were  computed  using 
a  layer  thickness  of  60  km,  av  erage  velocity  of  6  km  'sec  and  velocity 
standard  deviation  of  20%.  Ps  conversions  and  reverberations  start  to 
attain  amplitudes  seen  in  the  PAS  data. 

Ficure  11:  Envelope  stack  of  the  three-component  data  of  four  earthquakes  (see 
Table  1)  recorded  at  PAS.  The  envelope  is  shifted  10  seconds  for  display 
purposes.  Lines  snow  predicted  coda  levels  for  assumption  of  scattering  Q's 
of  100,  200,  300,  and  400  using  the  wholespace  energy  flux  model  (equation 
1 0b  Scattering  Q  at  PAS  is  approximately  200  to  300  but  the  coda  appears 
to  decay  slightly  faster  than  predicted  by  the  model. 

Figure  12:  Coda  envelope  stacks  for  the  ID  simulation.  The  "10%"  envelope  is 
the  stack  of  IS  vertical  and  radial  synthetic  seismograms  produced  by  9 
realizations  of  the  30  km  thick  layer  model  with  average  velocity  of  5.5 
km  sec  and  velocity  standard  deviation  of  10%.  The  "20%”  envelope  is  the 
stack  of  20  vertical  and  radial  synthetic  seismograms  produced  by  10 
realizations  of  the  60  km  thick  layer  model  with  average  velocity  of  6 
km/sec  and  velocity  standard  deviation  of  20%.  The  straight  lines  are 
least-squares  fits  of  the  coda.  Qs  values  are  those  inferred  from  the  zero 
lapse  time  intercept  of  the  linear  fits  (see  text).  Note  that  coda  decay  in 
the  ID  simulations  agree  with  the  assumption  that  coda  energy  follows  a 
diffusion  law  for  leaking  into  the  halfspace.  Coda  decay  is  controlled 
entirely  by  this  process  in  the  ID  simulations. 

Figure  15:  Three  component  data  recorded  at  SCP  (State  College,  PA)  on  the 

DWSSN  intermediate-period  system  (Table  4).  Long-period  noise  seen  on  some 
horizontal  waveforms  (e.g.,  08/12/85)  was  largely  removed  by  the  bandpass 
filter  used  in  this  study. 

Figure  14:  Envelope  stacks  of  the  PAS  and  SCP  data  showing  least-squares  fits 


of  a  line  through  coda  with  lapse  times  greater  than  10  seconds.  Qs 
inferred  front  each  lines’  zero  time  intercept  is  also  displayed.  Coda  excited 
at  PAS  attains  higher  levels  and  falls  off  more  slowly  than  coda  at  SCP. 

Figure  Al:  Amplitude  spectra  of  the  nominal  Benioff  1-90  response  (top)  and 
response  differences  (lower  curves)  assuming  30%  variation  in  damping  and 
free  periods  of  the  seismometer  and  galvonometer  of  the  system. 

Parameters  hs,  hg,  Ts,  and  Tg  are  the  seismometer  damping,  galvonometer 
damping,  seismometer  period  (sec)  and  galvonometer  period,  respectively.  See 
text  for  explanation. 

Figure  A2:  Comparison  of  Benioff  1-90  three  component  data  from  the  02/01/73 
event  (top)  and  the  12/28/73  event  (bottom).  Polarities  have  been  reversed 
for  the  vertical  (Z)  and  east-west  (E)  components  of  the  12/28/73  event  for 
comparison  purposes.  The  vertical  components  are  simple  showing  a  single 
impulsive  P  wave.  The  arrow  shows  the  location  of  the  major  Ps  conversion 
considered  in  the  this  study.  It  occurs  primarily  on  the  N  component  for 
the  02/01/73  event  and  on  the  E  component  for  the  12/28/73  event.  Likewise 
it  is  not  obvious  on  the  other  respective  horizontal  components  (right  side) 
showing  that  both  horizontal  instruments  respond  similarly  to  the  same  wave 
propagation  effect. 
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